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Abstract 

We study the asymptotic limiting function W{{G},q) = lim„^oo g)^/", where P{G,q) 
is the chromatic polynomial for a graph G with n vertices. We first discuss a subtlety in 
the definition of W{{G},q) resulting from the fact that at certain special points Qs, the 
following limits do not commute: lim^^oo hniq^g^ g)^^" 7^ limg^q^ lim„^oo g)^/". 
We then present exact calculations of W{{G},q) and determine the corresponding analytic 
structure in the complex q plane for a number of families of graphs {G}, including circuits, 
wheels, biwheels, bipyramids, and (cyclic and twisted) ladders. We study the zeros of the 
corresponding chromatic polynomials and prove a theorem that for certain families of graphs, 
all but a finite number of the zeros lie exactly on a unit circle, whose position depends on the 
family. Using the connection of P{G, q) with the zero-temperature Potts antiferromagnet, we 
derive a theorem concerning the maximal finite real point of non-analyticity in W{{G},q), 
denoted qc and apply this theorem to deduce that qdsq) = 3 and qdhc) = (3 + ^/E)/2 for 
the square and honeycomb lattices. Finally, numerical calculations of W{hc, q) and W{sq, q) 
are presented and compared with series expansions and bounds. 
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1 Introduction 



An important question in graph theory is the following: using q different colors, what is 
the number of ways P{G, q) in which one can color a graph G, having n vertices, such that 
no two adjacent vertices have the same color? The function P{G,q), first introduced by 
Birkhoff is a polynomial in q of order n and has been the subject of mathematical study 
for many years @-@]; reviews are Refs. @],citegraphs . Clearly, a general upper bound on a 
chromatic polynomial is P{G, q) < g", since the right-hand side is the number of ways that 
one can color the n- vertex graph G without any constraint. Consequently, it is of interest 
to study the limiting function 

Wi{G},q) = lim P(G,g)^/" (1.1) 

where the symbol {G} denotes the limit as n oo of the family of n-vertex graphs of type 
G. 

This limit has some characteristics in common with the thermodynamic limit in statistical 
mechanics, in which one defines a partition function at a given temperature T and external 
field H as Z = er^'^ (where the Hamiltonian Ti describes the interactions of the spins 
(jj, and (3 = {kBT)~^) and then, starting with a finite, usually regular, (i- dimensional n- 
vertex lattice G = A with some specified boundary conditions, one considers the reduced 
free energy (per site) / in the thermodynamic limit. 



e 



f = lim Z^/" (1.2) 



(Here, / is related to the actual free energy F hj f = ~I3F.) For spin models in the physical 
temperature range, < /3 < oo, the partition function Z is positive. In the case of the 
chromatic polynomial, for sufficiently large g, P(G, q) > 0. In both cases, one naturally 
chooses the real positive 1/ra'th roots in the respective eqs. ( |1.1D and (p..2[). 

Although the number of colors q is an integer in the initial mathematical definition of the 
chromatic polynomial, one may generalize g to a real or, indeed, complex variable. We shall 
consider this generalization here and study the function iy({G}, q) and the related zeros of 
P{G, q) in the complex q plane for various families of graphs G. For certain ranges of real 
g, P(G, q) can be negative, and, of course, when q is complex, so is P(G, q) in general. In 
these cases it may not be obvious, a priori, which of the n roots 

P(G,g)i/" = {|P(G,g)|i/V™>} , r = 0,l,...,n- 1 (1.3) 



to choose in eq. (|1 . 1| ) . Consider the function iy({G}, q) defined via eq. ( |1.1|) starting with q 



on the positive real axis where P(G, q) > 0, and consider the maximal region in the complex 
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q plane which can be reached by analytic continuation of this function. We denote this 
region as Ri. Clearly, the phase choice in ( p..3|) for g G -Ri is that given by r = 0, namely 
P(G,g)Vn = |P(G,g)|i/n. However, as we shall see via exactly solved cases, there can also be 
families of graphs {G} for which the analytic structure oiW{{G},q) includes other regions 
not analytically connected to and in these regions, there may not be any canonical choice 
of phase in (pT3|) . We shall discuss this further below. 

Besides being of interest in mathematics, chromatic polynomials P{G, q) and their asymp- 
totic limits W{{G},q) have a deep connection with statistical mechanics, specifically, the 
Potts antiferromagnet (AF) Denote the partition function for the (isotropic, nearest- 

neighbor, zero-field) g-state Potts model at a temperature T as Z = Y.{a„} ^^^'^ with the 
Hamiltonian 

n = -JY.^.^., (1.4) 

m 

where cTj = 1, g are Zg- valued variables on each site z G A. Define 

K = f3J , a = e^ (1.5) 

For the Potts antiferromagnet (J < 0), in the limit T — > 0, i.e., K — oo, the partition 
function only receives nonzero contributions from spin configurations in which at ^ cxj for 
nearest-neighbor vertices i and j, and hence, formally, 

Z{A,q,K=-oo)=P{A,q) (1.6) 

whence 

exp(/(A, q,K= -oo)) = W{A, q) (1.7) 

where G = A denotes the lattice. 

However, as we shall discuss in detail, the limit ( [I.ID , and hence the resultant function 
VrdG} = A, q) is not well-defined at certain special points without specifying further 
information. This constitutes a fundamental difference between the limits ( |1 . 1| ) and ( |1.2D ; in 
statistical mechanics, if the point Kq lies within the interior of a given physical phase, then 
f{A,K) is a (real) analytic function of K. Furthermore, in statistical mechanics the limit 
K Kq for a physical K, and the thermodynamic limit n oo (with the d-dimensional 
volume vold{A) oo) commute ||^ 

lim lim Z^/" = lim lim Z^/" (1.8) 

These limits still commute for complex K In contrast, the definition of W^(A, q) involves a 
further subtlety, since at certain special points g^ the following limits do not commute (for 
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any choice of r in eq. ( |1.3|) ): 



lim lim P(G,g)^/" ^ lim lim PiG.qf'' (1.9) 

n^oo q—>qs Q^Qs n-->oo 

As we shall discuss, the origin of this noncommutivity of limits is an abrupt change in the 
behavior of P{G, q) in the vicinity of such a point Qs] for q ^ Qs, P^G, q) grows exponentially 
as the number of vertices n in G goes to infinity: P(G, q) ~ a" for some nonzero a, whereas 
precisely at g = gs, it has a completely different type of behavior, which, in all of the cases 
considered here is P{G, g^) = Co(gs) where Co(g) may either be a constant, independent of n 
or may depend on n in a way that does not involve exponential growth, such as (—1)"- The 
set of special points {qg} includes g = 0, g = 1, and, on any graph G which contains at least 
one triangle, also g = 2; at these points, g^) = 0. It is also possible for P{G, g^) to be 
equal to a nonzero constant at g^. We shall discuss this further in Section 2. 

Before proceeding, we mention that, in addition to Refs. some relevant previous 

works are Refs. |[T^-[0. 

This paper is organized as follows. In Section 2 we discuss a subtlety in the definition 
of g), and in Section 3 we present some general results on the analytic structure of 

this function in the complex g plane. In Section 4 we give a number of exact solutions for 
W{{G},q) for various families of graphs {G} and calculate the resultant diagrams, showing 
the analytic structure of VrdG}, g) in the complex g plane. Section 5 contains a discussion 
of zeros of chromatic polynomials for various families of graphs and, in particular, a theorem 
on the location and density of zeros of g) for certain families G. In Section 6 we present 
a theorem specifying the maximal value of g, for a given lattice A, where the region boundary 
B for ly (A, g) crosses the real g axis, and we apply this to specific lattices. Section 7 gives 
numerical calculations of Vr(A, g) for the honeycomb and square lattices and a comparison 
with large-g series. Section 8 contains some concluding remarks. 



2 Definition of W{{G}, q) 

In order to discuss the subtlety in the definition of W{{G},q), we first recall the following 
general properties concerning the zeros of chromatic polynomials. First, for any graph G, 

P(G',g = 0)=0 (2.1) 

Second, for any graph G consisting of at least two vertices (with a bond connecting them), 

P(G,g=l)=0 (2.2) 
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Third, for any graph G containing at least one triangle, 



P(G',g = 2) = if GDA (2.3) 

These properties are obvious from the definition of P{G^ q), given that one must color adja- 
cent vertices with different colors. Since P{G, q) is a polynomial, each of these zeros at the 
respective values go = 0, 1 or 2 means that P{G, q) must factorize according to 

P(G,g) = (g-go)'(^")g(G,g) (2.4) 

where &(go) is a positive integer and Q{G, go) 7^ 0. One may distinguish two particular cases 
that occur for cases we have studied: (i) 6(go) = &o + bin] (ii) 6(go) = ^o, where 6o and h\ are 
integers independent of n. In the first, case. 



W{{G], q) = {q~ qot Q({G}, qf'^ (2.5) 



so that the two different orders of limits in ( |1.9| ) do commute. This type of behavior is 
observed for tree graphs, our first example below. However, in all of the other cases which 
we have studied, the second type of behavior (ii) holds. Hence for these families of graphs, 
as a consequence of the basic fact that 

lima;^/" = P ^ ^ ^ (2.6) 
1 if X = ^ ' 

and hence, lim„_,oo linig^gg (g — go)^(''°)/" = and lim^^gg lim„„,oo(Q' — go)^^'"'^''" = 1, the 
noncommutativity of the limits in eq. ( |1.9|) follows (for any value of r in (|1.3| )): 



lim lim P(G,g)^/" = (2.7) 

n— >oo q-^qo 

whereas 

lim lim P(G,g)i/" = lim lim Q(G,g)i/" ^ (2.8) 

q— >go »oo q^qo n^oo 

More generally, eq. ( |1 . 1|) is insufficient to define W{{G},q) not just in the vicinity of a 
zero of P(G, g), but also in the vicinity of any special point g^ where the asymptotic behavior 
of P{G, q) changes abruptly from 

P(G,g)~a" as n ^ oo (2.9) 

with a a nonzero constant, to 

P{G,qs) = const, as n ^ oo (2.10) 
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In case (ii) of eq. (|2.4| ), with b{qs) = Bq, one encounters this type of abrupt change in behavior 
with the constant in eq. (|2.10|) equal to zero. This is the origin of the noncommutativity of 



hmits in eq. ( |1.9|) at g = 0, 1 and, for G D A, at g = 2. However, this noncommutativity is 
more general and can also occur when the constant in eq. ( p.lOj ) is nonzero. An example is 
provided by the point = 3 on the triangular lattice; there are 3! = 6 ways of coloring a 
triangular lattice graph (with the technical provision that for finite triangular lattice graphs, 
one uses boundary conditions which do not introduce frustration). Denoting such a triangular 



lattice graph as tri„, it follows that P(trin, 3) = 6, which is of the form of eq. (|2.1CI| ) with a 



nonzero constant. For such cases, where the constant in eq. ( p.lO|) is nonzero, one has 



,^ v^,.//" 

n—>oo q—^Qs 

while 



lim lim |P(G,g)|i/" = 1 (2.11) 



lim lim |P(G',g)p/" = \a\ (2.12) 

q^qs n^oo 

where, in general, a ^ 1. Finally, the set of points {g^} also may include a continuous set 
comprising part of a region boundary, as will be discussed in theorem 1, part (e) below. 

Because of the noncommutativity ( |1.9|) , the formal definition ( |1 . 1|) is, in general, insuf- 
ficient to define W{{G},q) at the set of special points {g^}; at these points, one must also 
specify the order of the limits in (|1.9| ). One can maintain the analyticity of W{{G},q) at 
these special points g^ of P{G, q) by choosing the order of limits in the right-hand side of eq. 
(PD: 

W{{G},qs)n,^ ^ lim lim P(G,g)i/" (2.13) 

As indicated, we shall denote this definition as -Dgn? where the subscript indicates the order 
of the limits. Although this definition maintains the analyticity of W^dG}, g) at the special 
points g^, it produces a function iy({G}, g) whose values at the points g^ differ significantly 
from the values which one would get for P(G, g^)^/" with finite-n graphs G. The definition 
based on the opposite order of limits, 

W{{G},qs)D„, ^ Jim lim P(G,g)^/" (2.14) 

gives the expected results like W{G,qs) = for g^ = 0, 1, and, for G D A, g = 2, as well 
as W(trin,q = 3) = 1, but yields a function W{{G},q) with discontinuities at the set of 
points {gs}. In our results below, in order to avoid having to write special formulas for the 
points g^, we shall adopt the definition Dgn but at appropriate places will take note of the 
noncommutativity of limits (|1.9|) . 

As noted in the introduction, the noncommutativity of limits ( p..9|) and resultant subtlety 
in the definition of iy({G}, g) is fundamentally different from the behavior of the (otherwise 
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somewhat analogous) function e^^"^'^^ in statistical mechanics (where for this discussion, 
we consider a general statistical mechanical model and its reduced free energy, /, and do 
not restrict to the Potts model). The set {g^} includes certain discrete points lying within 
regions in the complex q plane where W{{G},q) is otherwise an analytic function. Now, 
considering the thermodynamic limit of a statistical mechanical model on a lattice G = A, 
one knows that for physical K, after the additive term {(/2)K + h is removed (where ( is 
the coordination number of the lattice A, h = jSH, and this removes the trivial isolated 
infinities in / at = oo and h = oo), i.e. after defining /(A, K) = {(/2)K + h + fr{A, K), 
the function fr{A,K) is analytic within the interior of a given phase. We should remark 
that in our studies of the properties of spin models generalized to complex-temperature, we 
have established that there may be singularities in thermodynamic quantities in the interiors 
of (complex-temperature extensions of physical) phases; specifically, we proved a theorem 



(theorem 6 in Ref. [|14|) that on a lattice with odd coordination number, the zero-field 
Ising model partition function vanishes, and the free energy / has a negatively divergent 
singularity, at the complex-temperature point z = —1, where z = e~^^^ and Kj = f3Jj is 
the Ising spin-spin coupling. For the honeycomb lattice, z = —1 lies on a phase boundary 
|T5| , but for the heteropolygonal lattice denoted 3 ■ 12^, 2; = — 1 lies in the interior of the 



complex-temperature extension of the ferromagnetic phase However, in the quantity 

analogous to W{A,q), namely, e^^^'^\ this singularity is a zero, not a discontinuity and 
furthermore it is not associated with the type of noncommutativity analogous to ( |1.9|) . The 
reason for this is that when one factorizes the (zero-field) partition function in a manner 
similar to eq. (^^ ), 

Z(A,2) = (2 + l)"/2^,(A,z) (2.15) 

(which defines Zs{A,z)), the exponent of the factor {z + 1) is proportional to n, as in case 
(i) for eq. ( p.4| ), so that the following limits commute: 

lim lim Z(A„,z) = lim lim Z(A„, ^) = (2.16) 



One should also contrast the noncommutativity (|1.9| ) with the very different type of non- 
commutativity which applies to a symmetry-breaking order parameter such as a (uniform 
or staggered) magnetization in a statistical mechanical spin model (above its lower critical 
dimensionality, so that it has a symmetry-breaking phase transition). Here, in both the 
symmetric, high-temperature phase and the low-temperature phase with spontaneously bro- 
ken symmetry, if one removes the external field before taking the thermodynamic limit, the 
magnetization vanishes: 

lim lim M(A, K,H) =0 (2.17) 
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whereas in the low-temperature, symmetry- broken phase {K > Kc), there is a nonzero 
magnetization in the thermodynamic hmit: 



hm hm M(A, K,H)j^O 



for K> Kr 



(2.18) 



However, this noncommutativity is quite different from that in eq. (|1.9| ): this is clear from 
the fact that, among other things, (|1.9|) can occur at a discrete, isolated set of special points 
Qs (as well as possibly a continuous set on a region boundary of type (e) in theorem 1 below), 
whereas the noncommutativity in eqs. (|2.17|) , (|2.18| ) occurs throughout the low-temperature, 
broken-symmetry phase of the spin model and, indeed, can be used to characterize this phase, 
with the spontaneous magnetization M{K, 0) (defined of course by the second ordering of 
limits, ( |2.18| )) constituting the order parameter. 



3 Analytic Structure of W{{G},q) 

As noted, we shall consider the variable q to be extended from the positive integers to 
the complex numbers. Although for a given graph G, P{G, q) is a polynomial and hence, a 
fortiori, is an analytic function of q, the function lydG}, q) which describes the n ^ oo limit 
of a given family of graphs {G} will, in general, fail to be analytic at certain points. These 
points may form a discrete or continuous set; if the set is continuous, it may separate certain 
regions of the complex q plane, which we denote Ri. We shall denote the boundary separating 
regions Ri and Rj as B{Ri, Rj) and the union of all components of regional boundaries as 
B = \Ji jB{Ri,Rj). On these boundaries, W{{G},q) is non-analytic. We shall illustrate 
this with exact results below. These regions are somewhat similar to complex-temperature 
extensions or complex-field extensions of physical phases in statistical mechanical models. 
However, there are also some fundamental differences. One of these is the noncommutativity 
of limits discussed in the previous section. Another is that in the case of ({(?}, q), it is not 
clear what would play the role of the physical concept of a order parameter characterizing 
a given phase (and its complex extension, such as to complex-temperature). Therefore, 
we shall use the terms "region" and "region diagram" rather than (complex extensions of) 
"phase" and "phase diagram". 

A question that arises when one considers the somewhat related regions (phases) in 
complex-temperature or complex-field variables for statistical mechanical spin models con- 
cerns the dimensionality of the locus of points where the reduced free energy / is non-analytic. 
Where the premise of the Yang-Lee theorem |I6| , p!7| holds (i.e., for physical temperature 
and Hamiltonians with ferromagnetic, but not necessarily nearest-neighbor, two-spin inter- 
actions, H = cTiJijO^j — HJ2i(^i on arbitrary graphs), it states that the zeros of Z in 
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the complex /i = e'"^^^ plane lies on the unit circle |/i| = 1 and hence, in the thermodynamic 
limit where these merge to form the continuous locus of points where / is non-analytic in the 
/i plane, this locus is one- dimensional. In the case of complex-temperature, taking the Ising 
model for illustration, the locus of points in the z = e~^^ plane is usually one- dimensional 
for isotropic spin-spin couplings J, but on the heteropolygonal 4 ■ 8^ lattice, it fills a two- 
dimensional area in this plane even for isotropic couplings [l^. In all of the exact results 



for W{{G},q) that we shall present below, the dimension of the continuous locus of points 
where W{{G},q) is non-analytic, is dim{i3} = 1. 
We next present a general theorem. 

Theorem 1 

Let G be a graph with n vertices and suppose that P{G, q) has the form 

Na 



P{G, q) = q{q - l)[c,{q) + c,(g)a,(g)"} (3.1) 



where Cj{q) are polynomials in q. Here Co(g) may contain n-dependent terms, such as (— 1)", 
but does not grow with n like a". This form and the additional factorization ( |3.4| ) are 
motivated by the exact solutions to be presented below. Note that the fact that P{G, q) is 
a polynomial guarantees that, for a given G, Na is finite. For a fixed q and, more generally, 
for a given region in the complex q plane, we define a term ai{q) to be leading if for q in 
this region |ci^((3')| > 1 and |a^(g)| > for all j ^ £. Without loss of generality, we can 

write ( |3.1| ) so that the aj{q) are different functions of q. Then our theorem states that (a) 
if Na > 2 and there exists some i such that |a^(g)| > 1 in a given region of the complex q 
plane, then if in this region, < 1, the term aj{q) does not contribute to the limiting 

function W{G,q); (b) if A^^ = 1 and Co(g) ^ 0, then if |ai| < 1, this term again does not 
contribute to W{G, q), which is then determined by co(g); (c) if A^a > 1 and ae{q) is a leading 
term in a given region of the q plane then (i) if this region is analytically connected to the 
positive real axis where P{G, g) > so that r = in (p..3|). 



W{{G},q) = aeiq) (3.2) 

while (ii) if (i) is not the case, then at least in terms of magnitudes, one has the result 

\W{{G},q)\ = \a,Xq)\ (3.3) 

(d) the regional boundaries B separating regions where different leading terms dominate are 
determined by the degeneracy in magnitude of these leading terms: |af(g)| = |a£/(g)|; (e) a 
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regional boundary can also occur where one crosses from a region where there is a leading 
term ai{q) to one where there is no leading term but there is a nonzero Co(g) ; this type of 
boundary is given by the equation |a.^(q')| = 1. 

Proof: (a) is clear since if < 1, then lim„^oo = 0, so in this limit it does not 

contribute to W{G, q), which is determined by the leading term, as specified in part (c). Part 
(b) follows by the same type of logic. Note that if there is no co{q) term and if iVa = 1 with 
|ai(g)| < 1, then in this case, ai still determines W{G,q); an example of this is provided by 
the region |g — 1| < 1 for the tree graphs T„ to be discussed below. Part (c) expresses the 
fact that in the limit n ^ oo, the contributions of subleading terms are negligible relative 
to that of the leading term, and hence the limiting function lydG}, q) depends only on this 
leading term. As one moves from a region with one dominant term a^(g) to a region in which 
a different term a^'(g) dominates, there is a non-analyticity in q) as it switches from 

W{{G},q) = a^{q) to W{{G},q) = ai/{q) for r = and similarly for nonzero r. This also 
proves (d). Statement (e) follows in a similar way. 

It is possible that P{G, q) contains no term of the form Cj{q)aj{q)^ but instead only the 
term co(g). Moreover, in the case where P{G,q) does contain such Cj{q)aj{q)^ terms, there 
may exist a region in the q plane where |aj(5')| < 1 for all j = 1, ...,Na. In both of these 
cases, W{{G},q) is determined by the remaining function co(g). 

If G contains one or more triangles A, then one may express P{G, q) in the form 



P{G,q) = q{q-l)iq-2){co{q) + Y.<^M^Mn} H G D A (3.4) 



Na 

E 

In this case, the same theorem applies, but with the further factorization ( ^.41 ) taken into 
account. 

The subtlety in the definition of lydG}, q) resulting from the noncommutativity (|1.9| ) is 
evident in the forms ( p.l| ) and ( ^.4| ). With the definition ( |2.13 ), if there is a leading term 
at go = 0, 1, and, for {G} ^ A, also = 2, then 

\W{{G},qo)\ = \a,{qo)\ (3.5) 

rather than zero, even though P{G, go) = at these points. If there is no leading term in the 
vicinity of a given go, i.e., if |aj(go)| < 1 for all j = 1, Na, then, if there is a Co(g) term, 

\W({G},q)\ = lim lim |co(g)p/" = 1 (3.6) 

5— >(jO n^oo 
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4 Exact Solutions for W{{G},q) 



In this section we calculate and discuss exact solutions for q) for various families {G} 

of graphs. We believe that these give some interesting insights into the analytic properties of 
such functions and also into the exact results obtained by Baxter for the triangular lattice. 
Unless otherwise cited, chromatic polynomials can be found, together with further properties 
of graphs, in Refs. 0-0. 

4.1 Tree Graphs 

A tree graph T„ is an n-vertex graph with no circuits and has the chromatic polynomial 
P(T„, q) = q{q — Using the procedure discussed in section 2, we choose r = in ( |1.3| ) 

and obtain 

Wi{T},q) = q-l (4.1.1) 

This applies for all q; i.e., W{{T},q) is analytic throughout the entire (finite) complex q 
plane. 

4.2 Complete Graphs 

An n-vertex graph is termed "complete" and denoted if each vertex is completely con- 
nected by bonds (edges) with all the other vertices. Thus, is the triangle, the tetra- 
hedron, and so forth. The chromatic polynomial is P{Kn, q) = 117=0 ~ -^^^ ^ given n, 
we may choose r = in ( [1.3|) by starting on the positive real q axis at a value q > n — 1. 
This yields 

Wi{K},q) = l (4.2.1) 

With our definition ( p.l3|) , W{{K},q) is analytic in the entire (finite) q plane. The zeros 
of P{Kn,q) are comprised by the set {go} = {0, 1, .■■,n — 1} and the noncommutativity of 
limits (|1.9|) occurs at each of these points. 

4.3 Cychc Graphs 

For the cyclic graph Cn, i.e., the n-circuit, the chromatic polynomial is 
P{Cn,q) = {q- l){(g - 1)""^ + (-!)"}• We find that the analytic structure of W{{C},q) 
differs in the two regions Ri and R2 consisting of q satisfying Ig — 1| > 1 and Ig — 1| < 1, 
respectively. The boundary B separating these regions is thus the unit circle centered at 
q = 1. These regions are shown in Fig. |I|. We calculate 

W{{C},q) = q-l for q e Ri (4.3.1) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Re(q) 



Figure 1: Diagram showing regional boundaries comprising B for 14^({C}, q) for cyclic graphs 
and zeros of P{Cn, q) for n = 19. 



For q G R2, the first term in curly brackets, (g — 1)" ^ — > as n — »• 00, and hence P{Cn, q) — > 
{q — 1)(— 1)"". This function does not have a limit as n 00. However, we can observe that 

\W{{C},q)\ = l for gei?2 (4.3.2) 

|Ty({C},g)| is, in general, discontinuous along B. For the choice r = in (|1.3|), it is 
continuous at g = 2 and has a discontinuity of 

\im Wi{C},q) - hm W{{C},q) = 2 (4.3.3) 

at g = 0. In Fig. |I| we have also plotted the zeros of P{Cn, q) for a typical value, n = 19. 
These will be discussed in Section 5. 
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4.4 Wheel Graphs 

The wheel graph {Wh)n is defined as an {n — l)-circuit C„_i with an additional vertex joined 
to all of the n — 1 vertices of C„_i (which can be thought of as the center of the wheel), and 
is naturally defined for n > 3. We find that the diagram describing the analytic structure of 
W{{Wh},q) consists of the two regions Ri and R2 defined respectively by |g — 2| > 1 and 
|g — 2| < 1, with the boundary B consisting of the unit circle 19 — 21 = 1. (These -Ri and R2 
should not be confused with the regions discussed in the previous subsection; we define the 
regions Rj differently for each family of graphs.) This diagram is thus similar to Fig. ^but 
with the circle moved one unit to the right; for brevity we do not show it. Using 

P{{WhU q) = q{q - 2){(g - 2)"-^ - (-1)^ (4.4.1) 

we calculate 

W{{Wh},q) = q-2 for q e Ri (4.4.2) 

For g G -R2, since the first term in curly brackets in eq. ( [4.4. 1|) , (g — 2)*^"^, goes to zero as n ^ 
00, one is left only with the second term, —(—1)"; this term does not have a smooth limit as 
n —>■ 00 and hence neither does the l/n'th power of this quantity. However, |H^({H^/i}, q)\ = 
1 for q G R2- Formally, one may choose the l/ra'th root such that W{{Wh},q) = —1 for 
q G R2- The noncommutativity of limits in eq. (|1.9|) occurs at the discrete points g = 0, 1, 2 
and, for even n, also at g = 3 since P{(Wh)n, n even, g = 3) = 0. More generally, for g 7^ 2, 
the n — s> 00 limit is not well defined on the circle |g — 2| = 1. 

One can also study wheel graphs with some spokes removed, which have been of recent 
interest |18|. Let us define the "cut" wheel {cWh)n/ as the n- vertex wheel graph with i 
consecutive spokes removed. For example, from an analysis of the specific case i = 2, we 
find the same boundary B, |g — 2| = 1, as for the asymptotic limit of the wheel graphs, and, 
furthermore, W{{cWh},q) = W{{Wh},q). 

4.5 Biwheel Graphs 

The biwheel graph f/„ is defined by adjoining a second vertex to all of the other vertices in 
the wheel graph {Wh)n-i and is naturally defined for n > 4. Here, we find that the diagram 
describing the analytic structure of W {{U} , q) consists of the two regions Ri and R2 defined 
respectively by |g — 3| > 1 and |g — 3| < 1 with B consisting of the circle |g — 3| = 1. The 
chromatic polynomial for these graphs is 

P{Un, q) = q{q - l)(g - 3){(g - 3)-^ + (-1)"} (4.5.1) 
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and from this we calculate 



W{{U},q) = q~3 for g G i?i (4.5.2) 

For q G R2, the term (g— 3)""^^ ^ as n — 00, and one is again left with a discontinuous term 
(— 1)". As before, one has in general, that for q G R2, \W{{U}, g)| = 1, and one can formally 
choose the l/ra'th root such that W{{U},q) = —1 in this region. The noncommutativity of 
limits in ( p..9|) occurs at the discrete points g = 0, 1, 2, 3, and also, if n is odd, at g = 4 since 
P{Un, n odd, g = 4) = 0. 

4.6 Bipyramid Graphs 

The bipyramid Bn is formed from the biwheel Un by removing the bond connecting the 
adjoined vertex to the center vertex of the biwheel. A bipyramid graph can be inscribed on 
the 2-sphere S"^ and, in this sense, can be considered to be 2-dimensional. The chromatic 
polynomial for i?„ is 

P{Bn, q) = - 2)"-^ + (g - l)(g - S)"-^ + {-If {q' - 3g + 1)} (4.6.1) 

Here we find a more complicated diagram describing the analytic structure (see also Ref. 
13); this is shown in Fig. |^ and consists of three regions: 

5 

Ri : Re{q) > - and |g - 2| > 1 (4.6.2) 



2 

5 
2 



5 

R2 : Re{q) < - and |g - 3| > 1 (4.6.3) 



and 

i?3:|g-2|<l and |g - 3| < 1 (4.6.4) 
The boundaries between these regions are thus the two circular arcs 

S(i?i, R,):q = 2 + e^\-^<9<^ (4.6.5) 

and 

27r Att 

B{R2, Rs):q = 3 + e'\—<<P<— (4.6.6) 
together with the semi-infinite vertical line segments 

BiR^,R2) = {q}: Reiq) = ^ and \Im{q)\ > ^ (4.6.7) 
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Figure 2: Diagram showing regional boundaries comprising B for W{{B},q) for bipyramid 
graphs and zeros of P{Bn, q) for n = 29. 

These meet at the intersection points g = 5/2 ± i\^/2. We find that 

W{{B},q) = q-2 for g G i?i (4.6.8) 

For the other regions, we have, in general, 

\W{{B},q)\ = = \q-3\ for q E R2 

= 1 for g G i?3 (4.6.9) 

With specific choices of 1/n'th roots, one can choose W{{B}, q) = q — 3 in R2 and —1 in R^. 

This example provides an illustration of the noncommutativity of limits (|1.9|) for a real 
non-integer point go, namely, go = (l/2)(3 + \/5) = Be^, = 2.618.., where the r'th Beraha 
number Bcr is given by [p!9| 

= 4cos2(7r/r) (4.6.10) 
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for r = 1, 2, ... The point Be^ lies in the region R3 where — 3g + 1 is the dominant term 
and is one of the two roots of this polynomial (the other root lies in region R2 and hence 
plays no role in W{{B}, q)). In Fig. |^ we have also plotted zeros of the bipyramid chromatic 
polynomial P{Bn, q) for a typical finite n = 29. These will be discussed in Section 5. 

Since Vr({G}, q) is bounded above by q, it is common to remove this factor and define a 
reduced function 

Wri{G},q)=q-'Wi{G},q) (4.6.11) 

which has a finite limit as |g| —>■ 00. There have been a number of calculations of Taylor 
series expansions in the variable l/(g — 1) for functions equivalent to Wr in the case where 
G is a regular lattice; see, for example, Ref. (and earlier references therein) for the 
square, triangular, and honeycomb lattices. Clearly, these series expansions rely on the 
property that, for these lattices, Wr{{G},q) is an analytic function in the 1/q plane at the 
origin, 1/q = 0. This analyticity of Wr{A,q) at 1/g = is proved by exact results for the 
triangular lattice and is strongly supported by numerical calculations of zeros of P(A, q) for 



A = sq,hc [|T2|, which show that the respective regional boundaries for these three lattices 
are compact, and do not extend to infinite distance from the origin in the complex q plane. 
However, our exact result for the infinite-n bipyramid function W{{B},q) and its region 
diagram demonstrates that, in general, the infinite-n reduced limit Wr{{G},q) of chromatic 
polynomials for a given family of graphs {G} is not guaranteed to be analytic at 1/g = 0: 
in the case of the bipyramid graphs, the portion of the regional boundary B comprised by 
the line segment ( [4.6.7| ) runs vertically through the origin of the 1/q plane, and Wr{{B}, q) 
is not analytic at 1/g = 0. 

4.7 Cyclic Ladder Graphs 

The cyclic ladder graphs with 2n vertices can be visualized as two n-circuit graphs (rings) 
Cni one above the other, with the z'th vertex of one n-circuit connected by a vertical bond 
to the i'th vertex of the other n-circuit. The chromatic polynomial for this family of graphs 
was calculated in Ref. |2ll] (where they are called prism graphs): 



P(L2„, q) = {q' -3q + 3)" + (g - 1){(3 - g)" + (1 - g)«} + - 3q + 1 (4.7.1) 

From this we compute the region diagram shown in Fig. ^ consisting of four regions: -Ri, 
R2, R*2, and R3 in which, respectively, (1) (g^ - 3g + 3)", (2)-(2)* (1 - g)", and (3) (3 - g)" 
are the leading terms. We find 

iy({L},g) =g2-3g + 3 for q E Ri (4.7.2) 
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\W{{L},q)\ = \l~q\ for q e R2 or R* (4.7.3) 

\W{{L},q)\ = \3-q\ for q e {4.7 A) 

Formally, one can choose 1/n'th roots so that q) = 1 — q and 3 — g in the respective 

regions (i?2,-R2) -^3- Note that, in contrast to the situation with the bipyramid graphs Bn, 
there is no region where all \aj{q) | < 1 so that the Co(g) term (which is equal to g^— 3g+l here) 
never dominates. The boundary between regions R2, R2, and R^ is the line |1 — g| = |3 — g| 
where these are both leading terms, which is comprised of the line segments 

B{R2, R3) = {q} : Re{q) = 2 , < Im{q) < ^2 (4.7.5) 

and its complex conjugate for B{R2,Rz)- Similarly, the boundary separating Ri from R^ is 
the locus of solutions to the degeneracy condition |(?^ — 3g + 3| = |3 — g| where these are 
leading terms. This boundary runs vertically through the origin g = and extends over, on 
the right, to two complex-conjugate triple points g = 2 ± \/2i where it meets the vertical 
line boundary ([4.7.5|) and its complex conjugate. Similarly, the boundaries separating Ri 
from R2 and Ri from i?2 are comprised of the locus of solutions to the degeneracy condition 
|g^ — 3g + 3| = |1 — g| where these are leading terms; as shown in Fig. ^ these boundaries 
extend from the above triple points over to a four-fold intersection point at g = 2. In passing. 



we note that our region diagram differs from that reported in Ref. [^, where the right- 
most curves were thought to terminate and hence not completely separate from Ri the two 
additional regions that we have identified as R2 and i?2- Fig. ^, we also show zeros of 
P{L2n, q) for 2n = 38. We shall discuss these in Section 5. 

4.8 Twisted Ladder (Mobius) Graphs 

One may also consider ladder graphs the ends of which are twisted once before being joined; 
these graphs are denoted twisted ladder or Mobius graphs, M2„. (It is easy to see that if one 
twists the ends an even number of times, this is equivalent to no twist, and any odd number 
of twists are equivalent to a single twist.) The chromatic polynomial is the same as that for 



L2n except for the co(g) term pi 



P(M2„, g) = (g^ - 3g + 3)" + (g - 1){(3 - g)" - (1 - g)"} - 1 (4.8.1) 

Since there is no region where the constant term co(g) (equal to —1 here) is dominant, we 
find that 

W{{M},q) = W{{L},q) (4.8.2) 
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5 Theorem for Zeros of Chromatic Polynomials for 
Certain {G} 

A general question that one may ask about zeros of chromatic polynomials is whether all, 
or some subset, of the zeros for an n- vertex graph G in the family {G} lie exactly on the 
boundary curves B. One knows that as n ^ oo, aside from the discrete general set of zeros 
of P{G,q), viz., go = 0,1, and, for graphs containing one or more triangles, q = 2, the 
remainder of the zeros merge to form the union of boundaries B separating various regions 
in the complex q plane. (Some of the set {qo} may also lie on B.) However, the fact that the 
zeros move toward, and merge to form, this boundary B in the n ^ oo limit does not imply 
that, for finite graphs G, some subset of zeros will lie precisely on B. We have investigated 
this question and have found that there do exist some families of graphs {G} for which the 
zeros of P{G, q) (aside from certain members of the set {qo}) he exactly on the respective 
boundary curves B. We shall present a theorem and proof on this. Interestingly, we find 
that in all such cases, B consists of a unit circle centered at a certain integral point on the 
positive real q axis. We emphasize, however, that this type of behavior is special and is not 
shared by other families of graphs that we have studied. Furthermore, for the families of 
graphs for which the theorem does hold, the positions of the unit circles differ for different 
{G}. Finally, as we shall show, the zeros populate the full circle with constant density. 

We find that the theorem applies for the following three families of graphs: (i) cyclic; 
(ii) wheel, and (iii) biwheel. We begin with the cyclic graphs. The form of P{Gn, q) differs 
depending on whether n is even, say n — 2m or odd, say n — 2m + 1. For even n > 4, we 
calculate the factorization 

' f (2j+l)7ri 1 r -(2i+l)7ri ~, 

P{Cn=2m. q) = q{q - 1) R " (1 + e^^)}{q - (1 + e^^)} (5.1) 

j=0 

and for odd n > 5, 

m— 1 . . 

P{Cr,=2n.+i, q) = q{q - l){q - 2) H {q - {1 + e^)}{q - (1 + e^)} (5.2) 

Special cases for lower n are P{C2, q) — q{q — 1) and P{Cz, q) — q{q — l){q — 2). 
For the wheel graphs, for odd n > 5, we find the factorization 

PiiWhU^m+i, q) = q{q -l){q- 2) - (2 + e^)}{q - (2 + e"^)} (5.3) 

and for even n > 6, 

m-2 

Pmh)n=2m. q) = q{q - l){q - 2){q -3)l[{q-{2 + e^)}{g - (2 + e^)} (5.4) 

3=1 
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Special cases for lower n are P{{Wh)'i, q) = q{q — — 2) and P{(Wh)4, q) = q{q — l)(g 
2)(g-3). 

For the biwheel graphs we calculate for even n > 6, 



P{Un=2m, q) = q{q -!)(?- 2){q - 3) {q - (3 + e^^)}{q - (3 + e"^^)} (5.5) 



j=0 

and for odd n > 7, 

m-2 



P{Un=2m+i, q) = q{q - l){q - 2){q - 3)(g - 4) J] - (3 + e^)}{g - (3 + e^)} (5.6) 



Special cases for lower n are P^U^, q) = q{q — l)(g — 2)(g — 3) and P{U^, q) = q{q — l)(g — 
2)(g-3)(g-4). 

These factorizations constitute a proof of the following 

Theorem 2 

Except for isolated zeros at g = 1 for C„, at g = 0, 2 for (Wh)n, and at g = 0, 1, 3 for Un, the 
zeros of P{Cn, q), P{{Wh)n, q), and P{Un, q) all lie on the respective unit circles |g — g©] = 1 
where g0(C„) = 1, qQiWhn) = 2, and qQ^Un) = 3. Furthermore, the zeros are equally spaced 
around the respective unit circles, and in the n ^ oo limit, the density g{{G},6) of zeros 
on the respective circles g = g© + e*^, — vr < ^ < vr, is a constant, independent of 6. If one 
normalizes g according to 

" g{{G},e)de=l (5.7) 



then 

gi{G},e) = l- for {G} = {C}, {Wh}, {U} (5.8) 
/vr 

For the cyclic ladder and twisted ladder graphs and M2n, we find the type of behavior 
which occurs with complex-temperature zeros of spin models: the zeros lie close to, but not, 
in general, precisely on, the asymptotic boundaries B. This is illustrated by the plots of 
zeros of P(L38, g) in Fig. ^ As one also finds in calculations of complex-temperature zeros 



in statistical mechanical spin models (see, e.g., p2|, ^]), the densities of zeros along certain 
boundary curves are very small; in Fig. ^ this occurs on S(i?i,i?2) near the intersection 
point g = 2. Similar low densities of zeros were observed for P{tri, g) on the boundary near 
g = and the right-most boundary near g = 4 |]T^. In statistical mechanics, the density of 
zeros g near a critical point Zc behaves as 

iz-z^r-"' (5.9) 
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where a' denotes the critical exponent describing the (leading) singularity in the specific heat 
at z = Zc as one approaches this point from within the broken-symmetry phase: Cgmg ~ 
\z — Zc\~"' p2|. Equivalently, the (leading) singularity in the free energy at z = Zc is given 



by fsing ~ 1-2 — Zc\'^~°'' . (Similar statements apply for the approach to z^ from within the 
symmetric phase with the replacement a' — > a.) Analogously, in the present context, the 
density of zeros of P{G, q) for a finite-n graph G near a singular point is determined by the 
nature of the singularity in the asymptotic function iy({G}, q): if one denotes the singularity 
in the function In VTdG}, q) at a point qc as 



In W{{G}, q)smg ~ |g - gcr"^ (5.10) 

where a'^ will, in general, depend on the direction (0) of approach to qc, then the corre- 
sponding density of zeros of P{{G}, q) as one approaches this point is 

(7({G},g)~|g-g/-"^ (5.11) 

For the bipyramid, as is evident in Fig. ^, one sees that, except for the general zeros at 
g = and 1 and a zero very near to g = Be^ = 2.618..., the inner zeros do lie near to the arcs 
forming the boundaries B{Ri, R3) and B{R2, R3), but the outer zeros do not lie very close to 
the line segments of B{Ri, R2), given by eq. ( [4.6.7| ), and only approach these line segments 
slowly as n increases. Since this latter behavior only occurs for the part of B extending to 
q = 5/2 ± ioo, it is plausible that it may be connected with the fact that this component 
of the boundary is noncompact. This inference is also consistent with the fact that of the 
families of graphs which we have studied, the bipyramid graphs form the only family with a 
noncompact B and the only family for which we have observed this deviation. 

6 Theorem on Singular Boundary of W{A, q) 



In recent work @-|25| on complex-temperature and Yang-Lee (complex-field) singularities 
of Ising models, it has been quite fruitful to carry out a full complexification of both the 
temperature-dependent Boltzmann weight u = z'^ = e~^^ and the field-dependent Boltz- 
mann weight /i = e~'^^^ and to study the singularities in the two-dimensional manifold 
depending on {z,fi) or {u,fi). This approach unifies the previously separate analyses of 
complex-temperature and Yang-Lee singularities; one sees that a given singular point {zc, fic) 
or {uc, fJ'c) in the manifold manifests itself as a singular point in the complex z ot u plane 
for a fixed /i and equivalently as a singular point in the complex /i plane for fixed z or u. 
We find this approach to be equally powerful here. Starting from the relation (|1.6|) be- 



tween the chromatic polynomial and the T = Potts antiferromagnet on a graph G, we 
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consider the two-dimensional complex manifold spanned by (a, q) (where a was defined 



in eq. ( |1.5|) ). For sufficiently large g, namely, q > where ( is the coordination number 



of the lattice G = A, the Dobrushin theorem implies ||26| that the Potts antiferromagnet is 



disordered, with exponential decay of correlation functions, at T = 0. As one decreases q, 
the AFM ordering tendency of the system increases, and, as q decreases through a critical 
value depending on the dimensionality d and lattice A, the model can become critical at 
T = 0, or equivalently, K = — oo. As one decreases q further, the AFM critical temperature 
increases from zero to positive values (i.e. Kc increases from — oo to a finite negative value). 
The critical value qc thus separates two regions in q: (i) the q > qc region, where the system 
is disordered at T = and (ii) an interval of q < qc where the system has AFM long-range 
order at T = (and for a finite interval < T < T^, where = J/{kBTc)). Now, using 
the relations ( [I.6|) , (|1.7|) and making the projection from the (a, g) space onto the real a 
axis, just as the disordered, Zjv-symmetric phase of the Potts AF must be separated by a 
non-analytic phase boundary from the broken-symmetry phase with AFM long-range order 
1 2^, so also, making the projection onto the real q axis for the W{A,q) function, it follows 



that the range q > qc and an adjacent interval q < qc must be separated by a non- analytic 
boundary. Furthermore, just as, by analytic continuation, the complex-temperature exten- 
sion of the disordered phase of the Potts antiferromagnet must be completely separated by a 
non-analytic phase boundary from the complex-temperature extension of the antiferromag- 
netically ordered phase, so also the region in the complex q plane containing the line segment 
q > qc must be completely separated from the region containing the adjacent interval to the 
left of qc in this plane. Since the zero-temperature criticality of the Potts AF and the critical 
value qc are both projections of the singular point {ac = 0, qc) in the manifold, we have 
derived the following theorem: 

Theorem 3: For a given lattice A, the point qc at which the right-most region boundary for 
W{A, q) crosses the real q axis corresponds to the value of q at which the critical point Oc 
of the Potts antiferromagnet on this lattice first passes through zero as one decreases q 
from large positive values. 

This point qc is the maximal finite real point of non-analyticity of Vr(A, q). 

We now discuss the application of this theorem to three specific 2D lattices. For this 
purpose, we recall that the g-state Potts model can be defined for non-integral as well as 
integral values of q because of the equivalent representation of the partition function 

Z=Y. v'^'''\<'''^ (6.1) 

G'CG 
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where G' denotes a subgraph of G = A, w = (a — 1), h{G') is the number of bonds and 
n{G') the number of connected components of G' . (Recall that one can see the connection 
of this with (|1.6|) by taking the K — — oo {v —1) limit of (|6.1|) , which yields the Whitney 
expression for P{G,q) 0.) 



6.1 Qc for the Honeycomb Lattice 

For the honeycomb (he) lattice, the paramagnetic-ferromagnetic (PM-FM) and PM-AFM 



critical points are both determined by the equation [|29| 



g^ + 3g(a- 1) - (a- 1)^ = (6.1.1) 



As q decreases in the range from 4 to g = Be^ = 2.618, one of the roots of eq. ( |6.1.1D 



increases from —1 to 0. This root can be identified as the AFM critical point a^iq) by going 
in the opposite direction, increasing q from its Ising value, q = 2 and tracking ac{q), which 
decreases from ac(2) = 2 — to a{q = Be^) = 0, where the AFM phase is squeezed out 
and there is no longer any finite-temperature AFM critical point, which now occurs only at 
T = 0. Hence, our theorem implies that 

q^(hc) = ^ = Be5 = 2.618.. (6.1.2) 

i.e., this is the value of q where the right-most regional boundary of W{hc,q) crosses the 
real axis in the complex q plane [0. We may compare this with the numerical calculation 



of zeros of P{hc,q), as a function of g, on finite honeycomb lattices in Ref. [12|. For this 
comparison, we first note that from the analogous study of the triangular lattice [0 , one sees 
the crossing point of a boundary curve increases by about Aq ~ 0.4 from an 8 x 8 triangular 
lattice with cylindrical boundary conditions (CBC's) to the thermodynamic limit. Assuming 
that a similar finite-size shift occurs for the honeycomb lattice, and noting that the curve of 
zeros calculated on the 8 x 8 he lattice with CBC's cross the real axis at g ~ 2.2, we find that 
these numerical results are consistent with the result of Theorem 3 for the thermodynamic 
limit. 

6.2 Qc for the Square Lattice 

For the square lattice, the PM-AFM critical point of the Potts antiferromagnet is given by 

(a + l)2 = 4-g (6.2.1) 
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a,(sg) = + (6.2.2) 
As q decreases from 4 to 3, this value of Oc increases from —1 to 0. Hence, 

qc{sq) = 3 (6.2.3) 

Our Theorem 3 then identifies qdsq) as the point where the right-most region boundary of 
W{sq, q) crosses the real axis in the complex q plane. Using the same rough estimate for the 
finite-size shift between the 8x8 square lattice with CBC's and the thermodynamic limit 
as was observed for the triangular lattice, viz., Aq ~ 0.4 and noting that a curve of zeros 



calculated for this finite square lattice crosses the real axis at g ~ 2.6 |12|, we see that our 



inference ( |6.2.3| ) is consistent with the numerical calculations in Ref. |T^. 



6.3 Qc for the Triangular Lattice 

Baxter's exact solution for W(tri, q) |12[ shows that in this case, 

q^{tri)=4 (6.3.1) 

where the right-most boundary, B{Ri,R2), crosses the real q axis. Theorem 3 implies that 
the two other singular points where the boundaries B{R2, R3) and Bi^R^, Ri) cross this axis, 
at g = 3.82... and q = 0, respectively, also correspond to singular points of the Potts 
antiferromagnet in the a plane at these two q values. 

7 Numerical Calculations of VF(A, q) for A = sg, he, tri 

The effect of ground state disorder and associated nonzero ground state entropy Sq has been 
a subject of longstanding interest. A physical example is ice, for which 5*0 = 0.82 ± 0.05 
cal/(K-mole), i.e., So/kB = 0.41 ± 0.03 [^, In statistical mechanical models, such a 
ground state entropy may occur in contexts such as the Ising (or equivalently, q = 2 Potts) 



antiferromagnet on the triangular |^ or kagome |3J] lattices, where there is frustration. 
However, ground state entropy can also occur in what is arguably a simpler context: one 
in which it is not accompanied by any frustration. On a given lattice, for sufficiently large 
q, Potts antiferromagnets exhibit ground state entropy without frustration; restricting to 
integral values of q, this is true for g > 3 on the square and honeycomb lattices, and for 



g > 4 on the triangular lattice For the given range of g on the respective lattices. 



since the internal energy U approaches its T = value U(T = 0) = —J{6aiaj)T=o = 
exponentially fast, it follows that for T ^ 0, i.e., K —00, limx^-oo PU = 0. Hence, from 
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the general relation 5" = (511 + /, it follows that, for this range of values of g, the ground 
state entropy (per site) and reduced free energy for the Potts antiferromagnet are related 
according to 

5o(A, q) = fiA,q,K = -oo) = In (A, q) (7.1) 

and hence lniy(A, q) is a measure of this ground state entropy. Accordingly, it is of interest 
to calculate Vr(A, q) for various lattices A and values of q. Moreover, from a mathematical 
point of view, for positive integer q, the numerical calculation of W^(A, q) gives an accurate 
measure of the asymptotic growth of P(A, q) ~ W^(A, g)" as the number of lattice sites 
n ^ oo. 

Extending our earlier calculation of W{hc, 3) we have calculated W{hc, q) for integer 
4 < g < 10. We use the relation for the entropy 

r/3 



S{I3) = S{(3 = 0) + PU{(3) - f U{(3')d(3' (7.2) 

Jo 



which is known to provide a very accurate method for calculating 5*0 [36|. We start the 
integration at /3 = with S{(3 = 0) = Ing for the g-state Potts antiferromagnet and utilize a 
Metropolis algorithm with periodic BC's for several L x L lattices with the length L varying 
over the values 4, 6, 8, 10, 12, 14, and 16 for all cases, and up to L = 24 for certain cases. Since 
U{K) very rapidly approaches its asymptotic value of as if decreases past about K = —5, 
the RHS of ( |7.2| ) rapidly approaches a constant in this region, enabling one to obtain the 
resultant value of S{P = oo) for each lattice size. For each value of q, we then perform a 
least squares fit to this data and extrapolate the result to the thermodynamic limit, and 
then obtain W from ( [7.1|) . We use double precision arithmetic for all of our computations. 
Typically, we ran several thousands sweeps through the lattice for thermalization before 
calculating averages. Each average was calculated using between 9,000 and 20,000 sweeps 
through the lattice. As we have discussed in Ref. for g > 3 on the honeycomb lattice 
(and also for g > 4 on square lattice), the finite-size dependence of 5*0 is not simply of the 
form 5o(A; L x L, g) = 5o(A, q) + c^j^\l-^ ; we fit our measurements with an empirical function 
of the form 

So{A; LxL,q) = 5o(A, g) + ci^^^L'^ + c%L-' + c^Il-' (7.3) 

As an example, we show in Fig. ^ the ground state entropy as a function of L~^, for the 
case g = 4 on the square and honeycomb lattices. As a check, we have confirmed that our 
measurements yield numbers consistent with the exact result 5*0 = for g = 2 and A = he, sq. 

In the course of these calculations, we have obtained a measurement of W{hc, 3) which is 
more accurate than, and in excellent agreement with, the value that we reported recently in 
Ref. (where we quoted the uncertainty very conservatively). This improvement is due to 
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Figure 4: Measurements of ground state entropy So, as a function of lattice size, for the 
q — A Potts AF on the honeycomb and square lattices. 
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q 


W{hc,q) 


he series 


W{sq,q) 


sq series 


3 


1.6600(5) 


1.6600 


1.53965(45) 


1.5396007.. 


4 


2.6038(7) 


2.6034 


2.3370(7) 


2.3361 


5 


3.5796(10) 


3.5795 


3.2510(10) 


3.2504 


6 


4.5654(15) 


4.5651 


4.2003(12) 


4.2001 


7 


5.5556(17) 


5.5553 


5.1669(15) 


5.1667 


8 


6.5479(20) 


6.5481 


6.1431(20) 


6.1429 


9 


7.5424(22) 


7.5426 


7.1254(22) 


7.1250 


10 


8.5386(25) 


8.5382 


8.1122(25) 


8.1111 



Table 1: Values of TV(A, q) for A = he, sq and 3 < g < 10 from Monte Carlo measurements, 
compared with large-g series. The entry for W{sq, 3) is from the exact expression. See text 
for further details. 



our use of larger lattices, up to 24 x 24, and double precision arithmetic, in the present work. 
Our results for W{he, q) are presented in Table |I| (with conservatively estimated uncertainties 
given in parentheses) and plotted in Fig. ^, together with the respective values obtained by 
evaluating the large-g series. For a lattice A, this series has the form 

W{A,q) = q(^f'W{A,q) (7.4) 
\ q / 

where, as above, ( is the lattice coordination number, and 

oo -1 

WiA,q) = l + Y.Wny\ y = (7.5) 

For the honeycomb lattice, W{he, qY = 1 + + 2y^^ + 4?/^^ + calculated through 0{y^^) 
20| . Because of the sign changes in the he series (the coefficients of the first five terms 



are positive, while those of the remaining four terms are negative), it is difficult to make a 
reliable extrapolation. Accordingly, for Table || and Fig. |^ we simply use a direct evaluation 
of the sum. As is evident from this figure, the agreement with our Monte Carlo calculation is 
excellent. From our result ( |6.1.2|) above, it follows that the large-g series cannot be applied 
below q = 2.62, and we have plotted it only down to the integer value, q = 3. For both the 
honeycomb and square lattices, we also show the exact results W{A, 2) = 1 for A = he, sq 
(as a superimposed circle and square) and W{A,0)£i^^ = W{A, l)D„g = (as a dot •), but 
we emphasize that these values assume the order of the limits in the definition Dnq in eq. 
( p.l4| ) and the respective values calculated with the other order of limits in the definition Dqn 
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in eq. (p.l3|) for these lattices (which values are not known exactly) could well be different 



from these Dnq values. 

It is also of interest to compare W^(A, q) for the other two regular 2D lattices, square and 
triangular. Unlike the honeycomb lattice, there have been previous Monte Carlo measure- 
ments for the square lattice for lattice sizes between 3x3 and 7x7 and q values up to 
10. We have extended these to considerably larger lattices, including 16 x 16 for all q values. 
As a check, for q = 3, from calculations on L x L lattices with periodic boundary conditions, 
for L = 4, 6, 8, 10, 12, 14 and 16, we obtain the fit 



So{sq, 3) = 0.431556 + 1.095289^-^ (7.6) 



This yields an asymptotic value which is in excellent agreement with the exact result |38 
So{sq,3) = (3/2)ln(4/3) = 0.43152311... (i.e., W{sq,3) = 1.5396007...); 



\So{sq, 3)exact " So{sq, 3)mc\ 



Soisq, 3) 



0.76 X 10"^ (7.7) 



exact 



The coefficient of the L ^ term agrees with a previous determination for this q = 3 case 
39| . Our fitting procedure for g > 4 has been discussed in conjunction with eq. (|7.3| ) above. 



We also compare our Monte Carlo calculations with the large-g series, ( [7.4|) , ( [7.5|) which, for 
A = sq was calculated to 0{y^^) in Ref. fS^]. The agreement is again excellent. In passing, 
we note that a calculation of the series for W{sq, q) to order 0{y^^) has been reported in Ref. 
P3|| , but we have checked that additional terms in this longer series have a negligible effect 
in the comparison of the series with our numerical results. Given our result that qdsq) = 3 
in ( |6.2.3|) , the large-g series cannot be applied for g < 3. It is interesting that the agreement 
between the series and our measurements is quite good even down to the respective region 
boundaries which we have deduced at qdhc) = 2.618 and qdsq) = 3. This suggests that 
the non-analyticities at these respective points on the honeycomb and square lattices are 
evidently not so strong as to cause the series to deviate strongly from the actual values of 
W{A,q). 

In passing, we remark that of course our results are consistent with the following rigorous 
bounds: (i) the general upper bound W{A, q) < q; (ii) the upper bound for the square lattice 
pO| , W{sq, q) < (l/2)(g — 2 + ^/(p^^^Aq+^) , which is more restrictive than (i) for q> 1/2; 
(ni) the lower bound applicable for any bipartite lattice, W{Aup.,q) > \/q — 1; and (iv) the 
lower bound for the square lattice W{sq,q) > {q^ — 3q + 3) / {q — 1). Note that for g = 2, 
both lower bounds (iii) and (iv) are realized as equalities, W{sq, 2) = 1. Although there is 
a range of q above 2 where (iv) lies below (iii), for g > 3, (iv) lies above (iii), i.e. is more 
restrictive. Some recent rigorous upper bounds on P{G, q) for general G have been given in 



27 



Ref. but these only improve the prefactor A multiplying P{G, q) < Aq"^, and hence still 
yield W{G,q) < q, as in (i). 

In the case of the triangular lattice, since, to our knowledge, there is no numerical evalu- 
ation in the literature of the expressions for W{tri, q) given in Ref. [1^, we have carried this 
out and plotted the resultant function in Fig. |^. The point g = 3 is an example of a special 
point qs discussed in section 2, where the behavior of P{G, q) changes abruptly from (|2.9| ) 
to ( p.lO|) and where, consequently, the two limits in (|1.9|) do not commute. As noted above, 
since there are just 6 ways of coloring a triangular lattice (equivalently, the ground state of the 
Potts AF is 6-fold degenerate), W(tri, 3)z)„^ = 1. We have indicated this with a symbol A in 
Fig. ^. However, with the other order of limits, (|2.13|) , lim^^s (trz, g) = W{tri,3)£)^^ ^ 1. 
(The actual value is iy(tri, S)^)^^ = 2.) In Fig. ||, at the other special points = 0, 1 (in- 
dicated with •) and gs = 2 (indicated with A) we have shown the values W{tr%,q)j:,^^ = 0. 
Again, however, because of the noncommutativity of limits in ( |1.9|) as discussed in section 
2, the value of W{tri,q)£)g^ calculated with the other order of limits, (|2.13|) is nonzero. At 
g = 1,2, it is positive, while at g = 0, the function has a discontinuity involving a flip in 
sign: 

lim W(tri, g) = — lim W(tri, q) (7.8) 

The sign of W{tri, g) for negative real g is unambiguous, since this interval is part of the 
region Ri and where consequently, there is a clear choice of 1/n'th root, given by r = 0, in 



( |1.3|) . This is also clear from the large-g series. In Refs. [|I2|, it was noted that the transfer 
matrix calculation used there can fail at the Beraha numbers g = Bcr] from the discussion 
that we have given in section 2 of this paper, we would view this as a specific realization of 
the general noncommutativity of limits ( |1.9| ). 

A general property one observes in Fig. |^ is that for these three lattices, for a fixed 
value of g in the range g > 4, W^(A, g) is a monotonically decreasing function of the lattice 
coordination number ( and a monotonically increasing function of the lattice "girth" 7, 
defined P| as the number of bonds, or equivalently, vertices contained in a minimum-distance 
circuit. Here, ( = 3,4,6 and 7 = 6,4,3 for A = hc,sq,tri. The dependence on the girth 
is easily understood: the smaller the girth, the more stringent is the constraint that no two 
colors on adjacent vertices can be the same. Concerning the dependence on (, we note that 
for tree graphs, where one can vary ( for fixed 7 (7 = 00), P{Tn, g) and hence iy({T}, g) are 
actually independent of ( (c.f. eq. ( |4.1.1|) ). This is also true for the magnitude |Vr(A, g)| for 
real negative g. 

For the bipartite (square and honeycomb) lattices, we observe the following general trend: 
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Figure 5: Plot of W{A, q) for A = sg, he, tri. At the special points = 0, 1 for all lattices, 
the zero values (denoted by the symbol •) apply for the order of limits in the definition Dnq 
in eq. ( |2.14| ). For the triangular lattice at the special points qs = 2,3, the respective values 
and 1 (denoted by the symbol A) also apply for this ordering of limits ( ^.14| ), as does the 
value W^(A, 2)0^^ = 1 for A = sq, he. The A = tri curve is plotted for the definition Dgn in 
eq. (CT ). 
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as q increases from g = 3, the ratios 



riy(A, q) 



iy(A,g) 



iy(A,g) 



(7.9) 



iy(A,g) 



max 



and the related 



rsi.K q) 



S{A,q;T = 0) 
S{A, q,T = oo) 



ln(iy(A,g)) 



(7.10) 



In q 



are monotonically increasing functions of q. The ratio rs{A, q) has a physical interpretation as 
measuring the residual disorder present in the g-state Potts antiferromagnet at T = 0, relative 



r5(/ic,3) = 0.461 and rs{sq,?,) = (3/2) ln(4/3)/ ln3 = 0.393 while rs{hc, 10) = 0.931 and 
rs{sq, 10) = 0.909. On the triangular lattice, as q increases from g = 4, one again finds that 

and monotonically increase; for example, rs{tri,A) = 0.273 while rs{tri, 10) = 0.852. 

Returning to the square and honeycomb lattices, although we cannot use our Monte 
Carlo method to evaluate W{hc, q) for negative g, we can use the large-g series, since the 
negative q axis is in the region Ri. Of course these series cannot be used all the way in to 
g = 0; in Fig. |^, we plot them up to g = —2. 

8 Conclusions 

In conclusion, we have presented some results on the analytic properties of the asymptotic 
limiting function IV({G},g) obtained from the chromatic polynomial P{G,q). We have 
pointed out that the formal equation ( |1 . 1| ) is not, in general, sufficient to define the function 
IVdG}, g) because of the noncommutativity of limits ( |1.9| ) at certain special points, and we 
have provided the necessary clarification for a complete definition of this asymptotic function. 
Using mathematical results on chromatic polynomials for several families of graphs {G}, we 
have calculated iy({G}, g) exactly for these families. From these results, we have determined 
the non-analytic boundaries separating various regions in the complex g plane for each of 
the VrdG}, g). We have also studied the zeros of chromatic polynomials for these families of 
graphs and have proved a theorem stating that for some families, all but a finite set of these 
zeros lie exactly on certain unit circles centered at positive integer points on the real g axis. 
Using the connection of chromatic polynomials to the partition function of the g-state Potts 
antiferromagnet on a lattice A at T = 0, in conjunction with a generalization to both complex 
g and complex temperature, we have presented another theorem specifying the position of 
the maximal (finite) real point gc(A) where M^dG} = A, g) is non-analytic and have applied 
this to determine gc on the square and honeycomb lattices. Finally, we have given Monte 



to its value at T = oo. This ratio is substantial; for exampL 



ie, from Table |l], one sees that 
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Carlo measurements of W{hc, q) (and W{sq, q)) for integral 3 < g < 10 and compared these 
with large-g series. Our results illustrate the fascinating and deep connections between the 
mathematics of chromatic polynomials and their limits on the one hand, and the statistical 
mechanics of antiferromagnetic Potts models on the other. 

This research was supported in part by the NSF grant PHY-93-09888. 
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